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^  LUMPED  GEOPOTENTIAL, HARMONICS  OF JORDER  29,  FROM  ANALYSIS 


OF  THE  ORBIT  OF  COSMOS  837  ROCKET 


/  (_J  H./ Hi  Her 


SUMMARY 


The  orbit  of  Cosmos  837  rocket  (1976-62E)  has  been  determined  at  36  epochs 
between  January  and  September  1978,  using  the  RAE  orbit  refinement  program  PROP6 
with  about  3000  observations.  The  inclination  was  62.7°  and  the  eccentricity 
0.039.  The  orbital  accuracy  achieved  was  between  30  m  and  150  m,  both  radial  and 
crosstrack. 

The  orbit  was  near  29:2  resonance  in  1978  (exact  resonance  occurred  on 
May  14)  and  the  values  of  orbital  inclination  obtained  have  been  analysed  to 
derive  lumped  29th-order  geopotential  harmonic  coefficients,  namely: 


9_0.2 
10  C29 


-  10  ±  15 


-  76  ±  12 


These  will  be  used  in  future,  when  enough  results  at  different  inclinations  have 
accumulated,  to  determine  individual  coefficients  of  order  29.  The  values  of 
lumped  harmonics  obtained  from  analysis  of  the  values  of  eccentricity  were  not 
well  defined,  because  of  the  high  correlations  between  them  and  the  errors  in 
removing  the  very  large  perturbation  (31  km)  due  to  odd  zonal  harmonics. 
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1  INTRODUCTION 

Cosmos  817  rocket,  I976-62E,  entered  the  following  orbit  on  1978  July  1 :  inclination 
82.75°,  perigee  height  440  km,  apogee  height  920  km,  period  98.4  min  and  eccentricity 
0.034.  Decay  is  expected  about  1984. 

The  orbit  seemed  promising  for  analysis  to  obtain  29th-order  'lumped'  geopotential 
harmonic  coefficients,  by  virtue  of  the  relatively  slow  passage  through  29:2  resonance 
with  the  Earth's  gravitational  field  -  wheti  the  satellite's  track  over  the  Earth's  surface 
repeats  every  two  days  while  the  satellite  makes  29  revolutions  of  the  Earth.  The  orbit 
has  been  analysed  for  the  nine-month  period  1978  January  to  September  when  the  effects  of 
29:2  resonance  were  significant:  exact  29:2  resonance  was  on  1978  May  14. 

The  orbit  was  determined  at  36  epochs  using  the  RAE  computer  program  RR0P6  ' . 

The  inclination  and  eccentricity  values,  after  clearance  of  perturbations,  were  fitted 

2 

by  least-squares  theoretical  curves,  using  the  THROE  program  .  The  best  fittings  gave 
lumped  29th-order  coefficients  which  are  nominally  the  most  accurate  so  far  obtained. 

These  lumped  coefficients,  with  others  at  different  inclinations,  will  be  used  to  obtain 
individual  29th-order  harmonic  coefficients.  A  few  other  values  have  previously  been 

3_5 

determined  ,  but  several  more  are  required,  especially  at  inclinations  near  90  . 

2  ORBIT  DETERMINATION 

2. 1  Observations 

Over  3100  observations  were  available  for  the  36  orbits  selected  for  determination, 
over  the  period  1978  January-September .  \bout  I 5Z  of  these  observations  were  rejected 
due  to  not  fitting  well,  to  leave  a  working  average  of  74  observations  per  orbit.  Three 
of  the  orbits  had  the  benefit  of  highly-accurate  Hewitt  camera  observations. 

The  largest  group  of  observations  used  was  about  2300  from  the  US  Navy;  a  further 
400  came  from  British  radar,  130  from  the  kinetheodoli te  at  the  South  African  Astronomical 
Observatory  (SAAO),  in  the  southern  hemisphere,  32  from  the  theodolite  in  Jokioinen, 
Finland  and  nearly  300  observations  were  from  volunteer  visual  observers,  supplied  by  the 
Anpleton  Laboratory  at  Slough. 

In  the  course  of  the  orbit  determinations,  it  was  noticed  that  several  groups  of  the 
South  African  observations  were  rejected,  many  with  consistent  errors  of  -3  or  -4  seconds. 
Correction  of  these  timing  errors  led  to  satisfactory  acceptance;  subsequently,  the  SAAO 
confirmed  that  such  errors  had  been  made  in  these  observations  (and  possibly  others)  over 
a  5-week  period  at  the  beginning  of  1978. 

2 . 2  The  orbits  and  their  accuracy 

The  36  computed  orbits  are  given  in  Table  I,  where  it  can  be  seen  that  the  sd  in 
inclination,  i  ,  varies  from  0.0003°  to  0.0017°,  the  rms  value  being  0.0008°;  by  compari¬ 
son,  the  rms  for  the  three  Hewitt  camera  runs  is  0.0005°.  For  eccentricity,  e  ,  the  sd 
—6  "“6 

varies  from  4  *  10  to  19  *  10  .  The  best  sd  values  of  i  and  e  are  equivalent  to 

30  m  in  position.  For  the  right  ascension  of  the  node,  .2  ,  the  average  sd  is  0.0009°, 
equivalent  to  100  m. 


Table  I 


Orbital  parameters  for  Cosmos  837  rocket,  with  standard  deviations 


[  ■ 

1 

|  KID 

Date 

1978 

a 

e 

i 

•• 

l 

43514.0 

Jan  6 

7056.6206 

0.038256 

62. 7406 

15.7015 

2.181 

2 

14 

11 

9 

10 

2 

522.0 

14 

7056.3017 

0.038355 

62.7417 

349.9948 

3.510 

2 

8 

8 

8 

8 

) 

532.0 

24 

7055.9888 

0.038501 

62.7420 

317.8549 

5.184 

3 

7 

8 

10 

10 

4* 

540.0 

Feb  1 

7055.6339 

0.038559 

62.7423 

292. 1450 

6 . 504 

3 

8 

3 

8 

10 

5 

548.0 

9 

7055.2032 

0.038676 

62.7386 

266.4234 

7.818 

3 

1  1 

1  1 

10 

10 

h 

555.0 

16 

7054.9588 

0.038812 

62.7382 

243.9120 

8.947 

3 

6 

6 

6 

7 

7 

564.0 

25 

7054.7220 

0.038966 

62.7377 

214.9628 

10.464 

3 

19 

1  7 

15 

17 

8 

574.0 

Mar  7 

7054.3304 

0.039081 

62.7388 

182. 7971 

12.042 

3 

13 

12 

10 

1  1 

9 

584.0 

17 

7053.8633 

0.039178 

62.74)1 

150.6240 

13.683 

2 

5 

8 

1  1 

8 

10 

594.0 

27 

7053.4390 

0.039300 

62.7375 

1 18.4402 

15.405 

3 

8 

13 

1  1 

15 

1 1 

604.0 

Apr  6 

7052.7617 

0.039338 

62.7345 

86.2450 

16.990 

4 

9 

8 

10 

10 

12* 

611  .0 

13 

7052.0397 

0.039385 

62.7365 

6 3  7034 

18.072 

4 

5 

5 

7 

8 

1  3 

616.0 

18 

705l..>292 

0.039357 

62. 7385 

47.5931 

18.934 

* 

1 1 

8 

9 

12 

14 

624.0 

26 

7050.8419 

0.039424 

62.7402 

21.8188 

20.228 

2 

1 ) 

9 

8 

1  1 

633.0 

May  5 

7050.1571 

0.019559 

62. 7371 

352.8086 

21.655 

6 

17 

14 

1  3 

15 

16 

641.0 

13 

7049.6910 

0.039619 

62.7373 

327.01 19 

22.885 

2 

7 

6 

5 

6 

17 

648.0 

20 

7049.36  31 

0.039675 

62.7356 

304.4)64 

24.008 

1 

8 

7 

6 

8 

15  i 

655.0 

27 

7049.0539 

0.039729 

62.7389 

281 .85b4 

25.127 

l 

4 

5 

5 

5 

14  1 

665.0 

Jun  6 

7048.5968 

0.039820 

62. 7417 

249.5988 

26.683 

1 

3 

1  3 

l  1 

9 

13 

JO  ( 

67  3.0 

14 

7048. 3420 

0. 039880 

62.7421 

223.7865 

27.95U 

1 

8 

6 

7 

8 

21 

bHl.O 

22 

7047.9572 

0.039976 

62.7381 

197.9700 

29.175 

5 

12 

8 

9 

14 

2  2 

688.0 

29 

7047.3403 

0.040027 

62.7362 

175.37)5 

30.235 

1 

9 

5 

4 

9 

2  3 

695.0 

Jul  6 

7046.7772 

0.0400)4 

62.7372 

152. 769  3 

31.344 

2 

15 

9 

8 

17 

2« 

701  .0 

12 

7046.3  310 

0.0400)4 

62. 7387 

133.3921 

32. 197 

3 

14 

6 

12 

16 

25  1 

707.0 

18 

7045.770) 

0.040048 

62. 7407 

1 14.0078 

33.181 

4 

16 

10 

16 

19 

( 

714.0 

25 

7045 . 1 294 

0.040104 

62.7)88 

91 . 385  3 

34. 349 

4 

1 1 

8 

9 

13 

27  | 

718.0 

29 

7044.8866 

0.0401  38 

62.7)84 

78.4571 

34.980 

2 

7 

5 

6 

9 

7  25  .0 

Aug  5 

7044.5810 

0.040210 

62.7  179 

55.8278  ' 

35.99J 

2 

5 

6 

<  | 

5 

7  )2.0 

12 

7044.  1127 

0.040245 

67.7)63 

3  3.1901  | 

37.117 

2 

9 

6 

9  i 

1  1 

7  )9.0 

19 

7044.1 | f 1 

0.04032) 

62.7377 

10.5516 

)8 . 192 

3 

12 

7 

10  i 

16 

746  .0 

26 

704  3.9807 

0.04041 1 

62.7402 

347.9142  ; 

)9 .285 

i 

3 

1  1 

8 

i  1  i 

16 

i 

754.0 

Sep  3 

7043.7319 

0.040486 

62.7414  i 

| 

40.489 

2 

8 

6 

0 

12 

761  .0 

10 

704 . .4281 

0,0406  J2 

62.7380 

299.396)  1 

41.587 

2 

9 

6 

7  j 

1  3 

u  1 

76  7  .0 

16 

7043.0906 

0.040578 

62.7)38  | 

279.98)5 

m2. 52b  i 

j 

2 

9 

6  { 

h 

!■  1 

77  3.0 

22 

7042.6663 

0.040587 

62.73J4  j 

260.5642  143.428 

! 

2 

8 

7  ! 

8 

i )  ! 

7  79.0 

28 

7041 .9341 

0.040579 

62.7344 

241.1  195 

44 , 386  | 

_ L 

1 

) 

10 

7  1 

7 

!  5  | 

14.424 

10 

81 .07) 
8 

257.305 

)0 

329.042 

10 

44.588 

10 

248.384 

7 

204.080 

16 

38. 375 
)  I 

237.579 

8 

81.595 

14 

291.491 

9 

1  56 . 1  30 
8 

268.728 
)  I 

22.  365 
10 

22.026 

1 5 

147.009 

6 

34.427 

7 

284.325 

5 

) 82 . 74 1 
12 

320.727 

8 

101 .303 

11 
3.401 

4 

270. 317 

16 

296.754 

15 

326.465 

18 

246.029 

12 

150.641 

9 

75.884 

5 

3.  319 
1  1 

292.658 

15 

22  3.25  3 
15 

42.709 
1  1 

336.842 
1  2 


52  7  1.0814 
3 

527  3.4  390 

5273.7898 

3 

5274.1879 

3 

5274.6707 

3 

5274.9447 

3 

5275.2103 

4 

5275.6499 
3 

5276.1 743 
3 

5276.6500 

3 

5277.4101 

4 

5278.2209 

5278.7942 

7 

5279.5664 

2 

5280. 3355 
7 

5280.8593 

3 

6281.2275 

1 

5281.5755 

2 

5282.0896 

4 

5282.  1761 
1 

5282. 8085 

5 

5283.5021 


0.02374  | 0.52 

10 

0 . 0  I  96  I  jo. 44 
8 

0.01666 
10 

0.03183 
16 

0.02684 
I  5 

0.01785 

14 

0.01677 

15 


0.02  384 
14 

0.02422 

6 

0.02807 
I  3 

0.05014 

21 

0.06162 

20 

0.05279 

28 

0.04268 
I  I 

0.04120 

26 

0.02927 


0.02429 
8 

0. 02 22 3 

7 

0.02080 
16 

0.01 796 
6 

0.036  6  7 
18 

0.05260 
9 

5284.1356  0.04248 

3  1  3 

5284.6377  0.04846 

4  21 

5285.2688  0.0558) 

5  21 

5285,4400  0.04316 

4  29 

5286.2633  0.02764 

3  18 

5286. b074  Q.02497 

I  I 

5286.909  3  |  0.0197) 

2  9 

5287.1  )65  |  0.01  104  0.M  6.6 

3  .  16 


0.56 
0.50 
0.44 
0.61 
0.64 
0.53 
0.60 
0.77 
0.54 
0.67 
0.49 
0.47 
0.  58 
0.53 
O.  )8 
0.  3a 

0.50 
0.36 
0.58 
0.  32 
0.63 
0.60 
0.64 

0.5) 

0.40 

0.49 

0.47 


5287.28)4 

3 

5287.56  18 


528  7  .  'J(j56 
2 


25.4  7  )  !  5288. 2H55  *  0.OJ.U5  0.^7 


Kt:y:  MJD  Mortified  Julian  Day 
a  semi  major  axis  (km) 
e  eccentricity 
i  inclination  (deg) 

**  fight  ascension  of  node  (deg) 

-  argument  of  perigee  (deg) 

*  orbits  with  Hewitt  camera  observat i ons 


1  I 

76.614 

12 

1  31  .47  3 
14 


5288.  7634 


4289.5685 

3 


0.01  I b7 
Ih 

0.01967 

9 

0.02966 
I  5 


0.6  3 
0.*7 
0.45 


l  i 

0,05089  '0.47 
21  1 

0.08715  |  0.55 
15  1 

- - 1 _ 


8.9 
H .  9 
8.6 
7 .6 

7.6 

6.9 
9.  3 

8.4 

9.1 

9.4 

5.9 

5. 1 

5.7 

8.9 
7.0 
7.0 
7.0 
7.  J 
7  .0 

7.7 

5.9 

5.9 

6.9 
5.6 
6  .  3 
4  .6 

4.9 

6  .  9 

8.0 


.7 
7  ) 
82 
76 
HI 
70 
70 
75 

8a  j 

i 

I 


6.2 


h  .0 
5.0 
5.6 


mean  anomaly  at  epoch  (deg) 

M|  mean  motion,  n  (deg /day) 

third  coefficient  in  polynomial  tor  mean  anomaly 
■  measure  of  fit 

0  time  coverage  of  observations  (days) 

N  number  of  observations  used 
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2 . 3  Motion  of  perigee 

Since  the  inclination  is  close  to  63°,  the  perigee  moves  very  slowly,  at  less  than 
0.2  deg/day,  and  the  argument  of  perigee  increases  from  2°  to  44°  in  the  265  days 

between  first  and  last  epochs,  as  shown  in  Fig  1.  The  odd  zonal  harmonics  in  the  geo¬ 
potential  have  a  great  effect  on  eccentricity  when  i  =  62.74°:  the  amplitude  of  the 
oscillation  is  nearly  50  km,  so  a  decrease  in  perigee  distance  of  over  30  km  is  to  be 

expected  as  w  increases  from  2°  to  44°.  In  fact,  a(l  -  e)  decreases  from  6787  km  on 

the  first  orbit  to  6756  km  on  the  last.  This  corresponds  to  a  decrease  in  perigee 

height  over  a  spherical  Earth  from  409  km  initially  to  378  km  at  the  end.  Fig  1.  Over 
an  oblate  Earth,  the  corresponding  values  are  409  km  and  386  km. 

2 . 4  Observational  accuracy 

The  residuals  of  the  observations  on  the  first  32  orbits  are  summarized  in  Table  2 
for  stations  with  5  or  more  observations  accepted.  The  residual  of  an  observation  is 
a  combination  of  the  observational  error  and  any  error  in  the  orbit,  and  the  value  given 
in  the  Table,  the  rms ,  produces  a  bias  towards  the  larger  values,  which  for  visual  obser¬ 
vations  relative  to  the  stars  are  usually  observations  made  in  poor  conditions  of  seeing. 
For  these  reasons  the  capability  of  visual  observers  in  good  conditions  is  usually  reckoned 
to  be  about  half  the  rms  residual.  For  the  US  Navy  station  29,  the  angular  residuals  are 
geocentric,  and  need  to  be  multiplied  by  a  factor  of  about  5  for  comparison  with  the  other 
( topocentri c)  observations.  All  observers  with  at  least  one  observation  accepted  have 
been  sent  copies  of  their  residuals. 

Table  2 

Residuals  for  observing  stations  with  more  than  5  observations  accepted 


Station 

N umb e r  of 
observations 
accepted 

Rms  residuals 

Range 

km 

Minutes  of 

arc 

RA 

Dec 

Total 

1  US  Navy 

187 

1  .7 

1  .9 

2.6 

2  US  Navy 

130 

2.  1 

2.2 

3.0 

3  US  Navy 

1  34 

2.0 

1  .9 

2.8 

4  US  Navy 

141 

1  .6 

2.0 

2.6 

5  US  Navy 

161 

1  .9 

1  .  7 

2.6 

6  US  Navy 

175 

1  .6 

1  .8 

2.4 

29  US  Navy 

625 

0.5 

0.  3* *• 

0.4* 

4  1  4X  Cape  Town 

28 

2.3 

2.3 

3.2 

2122x  Malvern  5 

1  3 

1  .6 

1  .2 

2  .  1 

2 1 2 5^  Street 

7 

2.9 

1 

3.5 

2 1 55x  Bahrein  2 

8 

2.8 

4 . 4 

5.2 

2265x  Farnham 

6 

2.0 

1  .  3 

2 . 4 

2303  Malvern  Hewitt  camera 

9 

0.02 

0.02 

0.0  1 

2414-^  Bournemouth 

59 

3.6 

1.5 

5.0 

2420x  Willowbrae 

28 

1  .*» 

2.  1 

2.8 

2577  Cape  kinetheodolite 

75 

0.8 

0.7 

1  .  1 

4I68J'  Vries 

1 1 

4.0 

1.8 

5.5 

6702x  Jokioinen 

2  1 

1  1.0 

1.6 

4 .  7 

*  Geocentric 

*•  Visual  stations 
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1  THE  29:2  RESONANCE 

i . I  Analysis  of  inclination,  i 

The  theory  for  29:2  resonance  is  detailed  in  Refs  3-5,  where  all  the  parameters 
used  here  are  defined.  The  first  term  of  the  inclination  equation  is 


di 

dt 


/R  ,3° 

ih-(f)  (29-2C°Si)?30,29,14 


_0,2 

sin  4  +  C^g  cos 


(1) 


where  4  =  2(uj  +  M)  +  29(11  -  v) 

is  the  resonance  angle,  v  being  the  sidereal  angle.  The  values  of  4  and  4  are  given 

in  Fig  2:  at  exact  resonance  4=0.  The  use  of  two  extra  terms  in  (i)  -  taking  (Y,q)  = 

(1,0),  (1,1)  and  (1,-1)  in  the  notation  of  Ref  4  -  gave  indeterminate  results,  very 

probably  because  the  (4  t  us)  terms  interfered  with  the  4  terms  as  a  result  of  the  very 

slow  variation  of  us  .  The  fitting  of  i  was  therefore  made  with  equation  (1),  in 

2 

integrated  form,  using  the  THROE  computer  program  . 

Before  being  fitted  by  THROE,  the  36  values  of  inclination  from  Table  1  were  first 
cleared  of:  zonal  harmonic  and  lunisolar  perturbations  (combined  maximum  value  being 
0.0070°),  using  the  PROD  program8  with  numerical  integration  at  daily  intervals;  tesseral 
harmonic  perturbations,  determined  by  PROP  (maximum  0.0016°);  and  atmospheric-rotation 
perturbations  (maximum  0.0018°)  determined  within  THROE,  using  an  atmospheric  rotation 
rate  A  of  1.0,  which  gave  a  better  fitting  than  the  other  alternative  tried,  A  =  0.9. 
Earth  and  ocean  tide  effects  were  not  taken  into  account,  and,  in  recognition  of  this, 
the  sd  of  one  value  of  inclination,  on  orbit  4,  was  degraded  from  0.0003°  to  0.0005°. 

A  density  scale  height  H  of  60  km  was  used,  appropriate  to  a  mean  height  of  430  km 
(0.75  H  above  mean  perigee  height). 

The  36  modified  inclination  values  were  then  fitted  by  the  integrated  form  of 
equation  (I),  using  THROE.  In  the  first  fitting  the  measure  of  fit  e  was  1.5,  so 
the  nine  worst-fitting  values  were  degraded  successively,  two  by  a  factor  of  3  and 
seven  by  a  factor  of  2.  Fig  3  shows  the  final  fitting,  and  the  values  with  their 
original  standard  deviations,  the  nine  degraded  values  being  marked  by  horizontal  bars. 

The  values  of  the  lumped  coefficients  finally  obtained  were: 

0,2  0  2 

io9C29  =  -  10  t  15  ,  io9S29  “  -  76  ±  12  ,  (2) 

with  •  =  0.88.  These  were  within  I  sd  of  the  values  obtained  on  the  first  fitting. 

Fig  )  shows  that  the  overall  effect  of  the  resonance  was  to  increase  the  inclination  by 
about  0.003°. 

-0,2 

lor  l97l-b2E  the  numerical  expression  for  in  terms  of  the  individual 

coefficients,  as  obtained  from  the  RAE  computer  program  PROF,  is 


So  ,2v 


1S<  32,29  * 


'•48C34,29  + 


°-82C36,29  -  °-70C38,29  '  °-66C40,29  + 


(3) 


vith  the  same  equation  for  S  ,  on  replacing  C  by  S  throughout. 


0  70 


7 


! . 7  Analysis  of  eccentricity,  e 

In  attempting  to  analyse  the  variations  in  eccentricity,  the  36  values  of  e  from 
Table  I  were  first  corrected  for  lunisolar  perturbations,  using  PROD,  and  then  the  values 


of  M.,  were  modified,  being  replaced  by 


:V,LeJ  -  'V 

J  ( t  ,  -  t  ) 
n+ 1  n 


,  wlie re  (M  )  is  the  value 
I  n 


of  M|  on  the  nth  orbit,  at  epoch  t^  .  This  allows  for  the  integrated  effect  of  drag 
between  successive  epochs'*,  as  required  by  THROE.  (This  correction  is  not  significant 
in  fitting  i  .)  In  the  THROE  runs  the  scale  height  II  has  to  be  taken  at  a  height 
I . 5H  above  perigee,  and  a  value  of  65  km  was  used. 

The  appropriate  equation  for  fitting  the  variation  in  e  due  to  resonance  is 


"  l6K29,29,14|C29  sin(I  "  v)  ~  S29  Cos('  "  Uj) 


u-'«3  . 


'29,29, 13)  29 


sin($  +  us)  -  S„ 


taking  (y,q)  =  (1,1)  and  (1,-1).  Accordingly,  the  values  of  e  ,  after  removal  of  air- 
drag  and  zonal -harmonic  perturbations  within  THROE,  were  fitted  by  THROE  using  equation  (4) 
in  integrated  form.  Unfortunately  the  values  obtained  for  the  lumped  coefficients 
(C,S)'^'  and  (C,S)~g’^  were  indeterminate,  probably  as  a  result  of  being  highly 
correlated  because  a  varies  so  slowly.  Also  the  value  of  e  was  high  (3.0),  and,  even 
after  five  values  of  e  had  been  degraded  in  accuracy,  thus  reducing  t  to  2.1,  the 
values  of  the  lumped  harmonics  were  still  all  less  than  twice  their  sd. 

The  first  possible  escape  route  from  the  'correlation  trap’  is  to  ask  whether  one 
pair  of  terms  in  (4)  is  small  and  can  be  ignored.  For  1976-62E,  the  PROF  computer  pro¬ 
gram  gives 


C29,29  f’'l4C3l,29 


10‘80C  33,29  '  4-MC35,29  "  5*08S7.29 


C29,29  ”  4-38C3l,29  +  4  ‘  *7<’  i  3 , 29 


8C33 ,2° 


”  ‘ 4  1  ('  37  ,  29  '  j 


and  similarly  for  S  .  Equations  (5)  indicate  that  the  most  probable  value  for  the 


-1,3/  1,1 


ratio  C , 


C,n^  is  about  1,  and  similarly  for  S  .  Since  the  multiplying  factors 


outside  the  two  curlv  brackets  in  equation  (4)  have  numerical  values  of  0.44  and  0.76 

- 1  ,  1 

respect i ve  Iv ,  the  most  probable  situation  is  that  the  C  terms  are  ot  about  the  same 

--1,3  . 

magnitude  as  the  C  terms  in  <.->),  so  that  neither  can  legitimately  be  ignored. 

Another  escape  route  might  be  found  by  taking  the  magnitudes  of  the  two  terms  as 
equal.  Since  thev  are  opposite  in  sign,  the  equation  (4)  would  reduce  to  the  form 


cos 


sin  1 )  sin  • 


where  K  is  constant.  So  if  ..  increased  from,  say. 


and  f i t  a 


'  10  !I0°,  it  would  be  possible  to  take  a  constant  value  of  sin  ~ 

i , ,q)  =  (1,0)  variation.  But  unfortunately  sin  increases  steadily  from  0  to  0.7 
and  cannot  be  taken  as  constant.  Thus  w  is  too  nearly  constant  to  allow  separation 
oi  the  effects  of  the  lumped  coefficients,  but  not  constant  enough  to  allow  useful 
s i mp 1 i f i cat  ions . 

Even  if  one  of  these  escape  routes  had  been  open,  however,  there  would  still  be 
another  obstacle  to  face:  the  amplitude  of  the  oscillation  in  perigee  height  is  very  large 
for  I476-62E  (about  46  km),  and  even  the  latest  set  of  odd  zonal  harmonics^  may  well  have 
uncertainties  of  up  to  2  km  at  this  inclination.  Since  the  change  due  to  resonance  is 
likely  co  be  less  than  1  km,  it  is  not  possible  to  remove  the  odd  zonal  harmonic  perturba¬ 
tion  with  adequate  accuracy. 


This  difficulty,  together  with  the  correlation  between  the  lumped  coefficients, 
prevents  a  satisfactory  solution.  But  the  values  of  eccentricity  still  need  to  be 
allotted  some  fitted  curve  for  future  comparison;  so  it  was  decided  to  adjust  the  values 


of  the  odd  zonal  harmonics  by  altering  J 7  so  that  (a)  the  value  of  e  was  near  minimal 
and  (b)  the  values  of  the  lumped  harmonics  were  consistent  with  the  maximum  credible 


values,  found  by  inserting 
cal  values  of  the  terms. 


(c,Sh,29  =  l0_5/:l2 


in  (5)  and  taking  the  sum  of  the  numeri- 


(The  individual  coefficients  are  generally  less  than  10  !'• 


see,  for  example,  Ref  8  -  and  they  will  rarely  all  be  of  the  same  sign.)  This  gives 

max.  (C,S)  ~  300  <  10  ^  and  max.  (C ,S)  2g  ’  ^  —  150  *  10  ^  .  The  best  solution,  which 

gave  values  within  1  sd  of  these  limits  and  had  a  reasonably  low  r.  ,  namely  t  =  1.6, 

-9  -9 

was  obtained  by  changing  from  the  standard  value  in  THROE,  -326  *  10  ,  to  -306  *  10  . 

The  resulting  values  of  lumped  harmonics  were: 


9-1,1  9-1,1 

UTC^  =  340  ±  750  10  S29 

9-- 1 , 3  o_- 1 , 3 

10  c29  =  -  320  ±  430  10  S,? 


(6) 


Tig  4  shows  the  values  of  e  ,  with  their  original  sd,  and  the  fitted  curve.  In  the 
;itting,  three  of  the  standard  deviations  were  relaxed  by  a  factor  of  3,  and  nine  by  a 
factor  of  2.  These  twelve  values  are  marked  by  horizontal  bars.  The  fitting  is  fairly 
satisfactory,  and  it  is  probable  that  determinate  values  for  the  lumped  coefficients 
would  have  emerged  if  they  had  not  been  so  strongly  correlated  and  of  similar  magnitude. 

Since  the  variation  of  perigee  height  with  time  (Fig  1)  is  near-linear,  another  way 
of  removing  any  error  in  the  calculated  amplitude  of  the  odd  zonal  harmonic  perturbation 
would  be  to  include  a  linear  term  in  the  fitting.  This  was  tried,  but  the  numerical 
.'efficient  of  the  linear  term  was  less  than  its  sd,  and  the  general  appearance  of  the 
fitting  was  not  appreciably  altered. 

Finally,  the  first  five  values  of  e  ,  which  seem  rather  low,  were  omitted,  and 

i  fitting  with  only  the  second  pair  of  lumped  coefficients  was  made.  The  resulting 

-9 

•  lines,  with  .1  =  “  296  *  10  ,  were: 


9 


9-“  » *  u_- 1 >  3 

10  C,9  =  -  4<S2  •  4  }  10  S29  =  202  4  46  *  (7) 

with  r  =  1.  34.  Tin*  standard  deviations  in  (73  are  mucli  lower  than  in  (6),  but  cannot  be 

-9 

accepted  as  realistic  because  a  change  ot  to  -2 76  *  10  changes  the  two  values  to 

-128  and  -49  respectively. 

The  most  useful  information  to  emerge  is  that,  in  all  the  21  fittings  attempted, 

_  ~  1 ,  3 

C  was  always  strongly  negative.  Since  its  numerical  value  seems  unlikely  to  exceed 

~  -9  -9 

150  *  10  ,  the  results  point  towards  a  value  somewhere  near  -150  x  10 

Although  the  present  analysis  of  e  has  not  yielded  good  values  of  lumped 
harmonics,  this  is  mainly  because  of  the  lack  of  knowledge  of  values  of  the  odd  zonal 
harmonics.  When  these  are  better  established,  a  successful  analysis  may  be  possible, 
since  the  orbital  data  are  basically  accurate  enough. 

4  COMPARISON  WITH  GEM  1 0B 

The  only  comprehensive  gravity  model  that  goes  to  order  and  degree  36  is  the 

9  .10. 

Goddard  Earth  Model  10B  and  recent  tests  with  accurate  resonant  orbits  indicate 

-9 

that  the  15th-order  terms  in  GEM  1 0B  are  probably  accurate  to  3  x  10  .  The  terms  of 

-9 

order  29  and  degree  30-36  are  likely  to  be  less  accurate,  and  an  error  of  5  *  10  may 
be  tentatively  assigned.  (The  use  of  GEM  10C,  which  is  the  same  as  GEM  I0B  to  order  and 
degree  36  but  goes  to  degree  180,  is  not  very  useful  here,  because  the  terms  of 
degree  >36  do  not  greatly  influence  the  lumped  coefficients.) 

-9 

Assuming  an  error  of  5  *  10  and  ignoring  the  terms  of  degree  38,40,...,  GEM  I 0B 

gives 

9-0,2  9-0,2 

10  C9g  =  -  8  ±  1 5  10  S  =  -  II  i  15  .  (8) 

Comparison  with  the  values  (2)  from  1976-62E  shows  good  agreement  for  C  and  disagreement 

-0,2 

for  S  .  The  high  negative  value  of  S?g  obtained  from  1976-62E  derives  directly  from 

the  main  increase  in  inclination  between  MJD  43615  and  43675  (see  Fig  3),  where  270°  < 

v  <  320°  (see  Fig  2)  and  sin  t  in  equation  (1)  therefore  has  a  value  between 

-0.64  and  -1.0.  This  'main  increase'  in  inclination  seems  securely  based,  because  it 

-0,2 

is  so  strongly  confirmed  by  all  the  outlying  points  in  Fig  3.  So  the  value  of 

from  equations  (2)  seems  preferable  to  that  from  equations  (8). 

_-l,3 

The  value  of  C?g  from  GEM  I0B  is  not  likely  to  be  very  realistic,  because 

of  the  neglect  of  terms  of  order  37,39,...;  but  the  7  =  33  term  in  GEM  10B  gives 

9--  I  ,  3 

a  contribution  of  -56  *  20  to  10  ,  which  is  consistent  with  the  strongly  negative 

value  indicated  by  I976-62E. 

5  CONCLUSIONS 

The  orbit  of  Cosmos  837  rocket  has  been  determined  at  36  epochs  spread  throughout 
the  first  nine  months  of  1978,  when  the  effects  of  29:2  resonance  were  being  felt.  About 
3100  observations  were  used,  of  which  about  757  were  supplied  by  the  US  Navy. 


10 


The  3b  orbits  obtained  are  given  in  Table  I  and  show  chat  the  sd  in  inclination 
varies  between  0.000  3°  and  0.0017°  while  the  sd  i"  "ccentricity  varies  irom  4  -  1 0  h  to 
l1'  *  10  The  best  standard  deviations  are  each  equivalent  to  30  ra  in  position. 

The  3b  values  of  inclination  (cleared  of  zonal-harmonic,  J  0  ,  air-drag  and 
lunisolar  perturbations)  were  fitted  with  a  least-squares  theoretical  curve  to  give  the 
following  values  of  lumped  29th-order  coefficients: 

9-0,2  9-0,2 
10  C.,9  =  -  10  *  15  10  S  =  -  7b  •  12  . 

The  analysis  of  eccentricity  was  not  so  successful  because  of  (a)  the  difficulty 

of  accurately  removing  the  very  large  perturbation  (31  km)  due  to  odd  zonal  harmonics, 

and  (b)  the  interference  between  the  two  pairs  of  lumped  harmonics,  caused  by  the  very 

slow  variation  of  the  argument  of  perigee  at  this  inclination  (62.74°).  However,  there 

— — 1,3  ^  —9 

is  a  strong  indication  that  C^g  is  strongly  negative,  of  order  -150  "  10 

The  values  of  lumped  harmonics  from  1976-62E  -  which  appear  to  be  the  best  so  far 
obtained  at  29:2  resonance  -  will  be  used  in  future  determinations  of  individual 
coef f icient s ,  when  results  for  a  variety  of  inclinations  are  available. 
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=  -10+15  and  10  S 
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=  -76+12 


These  will  be  used  in  future ,  when  enough  results  at  different  inclinations  have 
accumulated,  to  determine  individual  coefficients  of  order  29.  The  values  of 
lumped  harmonics  obtained  from  analysis  of  the  values  of  eccentricity  were  not 
well  defined,  because  of  the  high  correlations  between  them  and  the  errors  in 
removing  the  very  large  perturbation  (31  km)  due  to  odd  zonal  harmonics. 
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